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To finally resolve the controversial issue of whether or not the electronic structure of YbBe is 
nontrivially topological, we have made a combined study using angle-resolved photoemission spec¬ 
troscopy (ARPES) of the non-polar (110) surface and density functional theory (DFT). The flat-band 
conditions of the (110) ARPES avoid the strong band bending effects of the polar (001) surface and 
definitively show that YbBe has a topologically trivial B 2p-Yb 5d semiconductor band gap of ~ 0.3 
eV. Accurate determination of the low energy band topology in DET requires the use of a modified 
Becke-Johnson exchange potential incorporating spin-orbit coupling and an on-site Yb 4/ Coulomb 
interaction U as large as 7 eV. The DFT result, confirmed by a more precise GW band calculation, 
is similarly that of small gap non-Kondo non-topological semiconductor. Additionally the pressure- 
dependent electronic structure of YbBe is investigated theoretically and found to transform into a 
p-d overlap semimetal with small Yb mixed valency. 

PACS numbers: 71.20.Eh,71.27.+a,75.30.Mb 


A great deal of recent attention has been paid to the 
topological nature of strongly correlated systems, which 
include the topological Mott insulator [1, 2], the frac¬ 
tional topological insulator [3, 4], and the topological 
Kondo insulator (TKI) [5]. In these systems, the in¬ 
terplay between topological characteristics and strong 
electron correlations provides new interesting phenomena 
that can possibly be utilized for spintronic and quantum 
computing applications. 

The first candidate material for a TKI is SmBe, which 
has been predicted first theoretically [5-8], and then 
studied intensively by transport [9-11], angle-resolved 
photoemission spectroscopy (ARPES) [12-15], and scan¬ 
ning tunneling microscopy/spectroscopy (STM/STS) 
[16, 17] experiments to explore its surface states. Subse¬ 
quently, other 4/-electron systems have been proposed as 
TKTs and topological Kondo semimetals [18-24]. Two 
essential common ingredients for a non-trivial topological 
character are (i) band inversion between opposite parity 
4/ and bd states, caused by rare-earth mixed-valence, 
and (ii) a large spin-orbit coupling (SOC) provided by 
the 4/ states. At the simplest level a strongly correlated 
bulk topological insulator (TI) would have the generic TI 
property of protected, symmetry-required, spin-textured 
metallic Dirac cone surface states that span the insulat¬ 
ing bulk gap. 

YbBe of our present interest was proposed to be a 
TKI with the mixed-valence state of Yb being 2.2 (uf = 
13.8) based on the inverted Yb Af-bd bands obtained in 
the density-functional theory (DFT) + Gutzwiller band 
method [22]. However, early photoemission [25] and re¬ 


cent ARPES [26-29] show that the binding energy (BE) 
of the Yb 4 / 7/2 band is about 1 eV, indicating that there 
would be no f-d band inversion and so YbBe would not 
be a TKI. Then, inspired by the observation of (001) 
surface states having the appearance of Dirac cones [26- 
28], two ARPES groups proposed that YbBe would be a 
weakly correlated TI with band inversion between oppo¬ 
site parity Yb bd and B 2p bands [27, 28]. The topolog¬ 
ical origin of the observed surface states was questioned 
[29], however, because they were observed to not follow 
the expected linear Dirac cone dispersion and to exhibit 
time-dependent changes. Instead band bending and sur¬ 
face quantum well confinement arising from the (001) po¬ 
lar surface was suggested, while not explicitly proposing 
that YbBe is not a TI. 

The p-d band inversion TI scenario was supported the¬ 
oretically with DFT + SOC U {U = 4 eV) calculations 
[27, 30], but also with an incorrect 0.3 eV BE of the 
Yb 4 / 7/2 state and in contradiction to an earlier calcu¬ 
lation [31] using U = 7 eV that obtained a p-d inverted 
semimetal with the correct experimental Yb 4/ energy. 
These current experimental and theoretical uncertainties 
have prevented a consensus on the topological nature of 
YbBe. 

In this Letter, we report new ARPES experiments 
that definitively demonstrate the non-Kondo non-TI elec¬ 
tronic structure of YbBe and new DFT theory that agrees 
well with the experimental results and strongly supports 
the same conclusion. ARPES for the non-polar (110) 
surface reveals a clear p-d semiconductor gap with no 
in-gap surface states, whereas all surfaces of a TI sys- 
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tern must have surface states. Calculations incorporating 
the SOC and U into the modified Becke-Johnson (mBJ) 
potential [32] describe properly the BE of the Yb 4 / 7/2 
band and the observed ARPES spectra of a topologically 
trivial Yb Sd-B 2p band gap. We have also investigated 
the pressure-dependent electronic structure of YbBe and 
found that the high pressure phase is a topologically non¬ 
trivial p-d overlap semimetal with an Yb 4 / 7/2 BE of 
^0.5 eV, rather than an full insulator. This result ex¬ 
plains a recent experimental study of transport and Yb 
valence under pressure [33]. 

ARPES measurements were performed at the MER¬ 
LIN Beamline 4.0.3 at the Advanced Light Source in the 
photon energy {hv) range of 30—150 eV. An elliptically 
polarized undulator was employed, which allows selec¬ 
tion of s- and p-polarization of the incident light. A Sci- 
enta R8000 hemispherical electron energy analyzer was 
used with energy resolution set to « 20 meV [34]. Mea¬ 
surements were performed in a vacuum of better than 
5 X 10“^^ Torr for the sample cooled down to ^ 30 K. 

The band calculations were performed using the full- 
potential linearized augmented plane-wave (ELAPW) 
band method, as implemented in the WIEN2K package 
[47]. Eor the DET calculations, the PBE (Perdew-Burke- 
Ernzerhof) exchange-correlation functional was used in 
the GGA (generalized-gradient approximation). In the 
GGA -1- SOG + U method, a correlation energy oiU = 1 
eV was chosen to obtain the correct experimental value 
of the Yb 4/ BE of ~ I eV, which is consistent with 
the previous calculations [31, 48]. The mBJ potential 
is adopted to provide band gap corrections in agree¬ 
ment with the improved many-body but much more 
computation-demanding GW calculation [32, 49]. The 
details of the calculational methods are described in the 
Supplement [34]. 

Eor insulating hexaborides, the polarity of (001) sur¬ 
faces with different charge terminations can lead to n- 
and p-type band bending and quantum well states that 
make it difficult for ARPES to directly observe the bulk 
band gap. Also spectra from spatially inhomogeneous re¬ 
gions (z.e., both n- and p-type) can falsely appear to show 
Dirac cones or p-d overlap [29, 50]. While surface mod¬ 
ification and aging provide some control over the band 
bending and assist in the ARPES interpretations [50], 
these problematic band-bending effects can be avoided 
by instead measuring a non-polar surface such as the 
(110) surface whose charge neutral bulk-termination is 
schematically illustrated in Eig. 1(a). Eor this purpose, 
a (110) surface of YbBe was prepared from the natu¬ 
ral facet of a single crystal grown by the aluminum-flux 
method. After etching in hydrochloric acid and ion sput¬ 
tering of the surface, the sample was annealed to 1300 
°G in ultra high vacuum to produce a spatially uniform 
1x1 ordered surface [34]. 

X-point ARPES spectra measured along M-X-M at hu 
= 120 eV using two different linear polarizations of the 


incident light is shown in Eigs. 1(c) and (d). Above 
a strong Yb 4/ peak at —1.05 eV, the p-polarization 
spectrum shows a weak hole-band dispersion and a small 
electron-like intensity at Ep. A strong polarization se¬ 
lectivity of these states is revealed by the 5-polarization 
spectrum in Eig. 1(d), where the electron conduction 
state is totally suppressed and the valence hole band is 
strongly enhanced to manifest a triangular-like dispersion 
with a rounded-maximum and hybridization interaction 
with the Yb 4/ states. The strong hole-band intensity 
allows a quantitative fit (dashed line) to a two-band k • p 
non-parabolic dispersion model [34] with a band maxi¬ 
mum of 0.35 eV below Ep. 

Eigure 1(e) shows an enlarged view of the p- 
polarization spectrum in which the Yb 4/ spectral in¬ 
tensity tail has been divided out to obtain an enhanced 
image of the ^ 0.3 eV semiconductor band gap between 
the B 2p valence and Yb 5d conduction states. To further 
characterize the conduction band dispersion and energy 
minimum, K-dosing of the surface was used to induce a 
small n-type band bending until the electron pocket was 
increased in depth to 0.2 eV revealing enough of a dis¬ 
persion [34] to perform similar non-parabolic dispersion 
analysis. The resulting process exhibited no discernible 
surface band gap narrowing, thus allowing evaluation of 
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FIG. 1: (Color Online) (a) Schematic structure of the non¬ 
polar YbBe (110) surface, (b) Cubic BZ with (110) orienta¬ 
tion illustrating the locations of bulk X-points. (c,d) X-point 
spectra measured at hu = 120 eV with p- and s-polarization, 
illustrating the opposite polarization dependence of p-hole 
and d-electron states (e) Zoom of the p-polarization spectrum 
with the Yb 4/ spectral intensity removed to enhance the view 
of the ~0.3 eV band gap. Dashed lines are non-parabolic fits 
to the spectral intensity maxima (see text). 
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FIG. 2: (Color Online) Off-normal photon-dependent map of 
YbBe (110) at fixed kx = —0.53 A as shown in Fig. 1(b), 
of the spectral intensity of (a) valence band at —0.3 eV and 
(b) conduction band at Ef- The 3D bulk-like kz dependences 
of both the valence and conduction bands confirms the flat 
band conditions of the non-polar (110) surface. 


a band gap of 0.32 eV. 

To explicitly confirm that the bands shown in Fig. 1 
are bulk, we have measured their kz dependences in the 
process of locating the bulk X-points. Figures 2(a) and 
(b) show the ky-kz maps at fixed kx = —0.53A for the 
valence band at — 0.3 eV and the conduction band at 
respectively. Both constant energy cuts exhibit strong in¬ 
tensity features close to bulk X-points ai ky = 0 for hu 
= 76 eV and 120 eV as well as intensities at X-points of 
the second Brillouin zone (BZ) at ky = ±1.5A . The 
small vertical /i:; 2 -elongation of the X-point intensities in 
Fig. 2 is well accounted for by the inherent bulk band 
structure anisotropy (see Fig. 3) and the -broadening 
effect resulting from the finite inelastic mean free path of 
the photoelectrons. The pinning of Ep at the bottom of 
the conduction band is consistent with the negative sign 
of the bulk Hall coefficient [33, 51, 52], and consistent 
with flat-band conditions of the non-polar (110) surface. 
Hence both the valence and conduction bands shown in 
Fig. 1 are 3D-like bulk bands and do not originate from 
the 2D-like surface states. The strong polarization de¬ 
pendence in Fig. 1(d) also independently confirms that 
these states are not linear Dirac cone dispersions, which 
would instead exhibit some continuity of the same or¬ 
bital characters between the upper and lower parts of 
the Dirac cone. 

The bulk X-point spectrum in Fig. 1(e) exhibiting 
a clear small direct semiconductor gap between valence 
and conduction band states and the absence of in-gap 
surface states is the central experimental result of this 
study. The (110) ARPES definitively proves the absence 
of a p-d overlapping band structure and hence a lack of 
parity inversion that is the key first requirement for a 
topological electronic structure interpretation of previ¬ 
ous ARPES for the (001) surface. Therefore, the ob¬ 
served chirality in 2D surface states of YbBe (001) in 
circular-dichroism (CD) [26] and spin-resolved ARPES 
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FIG. 3: (Golor Online) DFT-Wien2k band structures of 
YbBe. (a) GGA + SOG ± U (7 eV) calculation yields a 
semimetallic p-d overlap with anti-crossing gaps, (b) mBJ + 
SOG ± U (7 eV) bands overlaid with open-core GW band 
results (dots). Both exhibit semiconductor band gaps. 


[28], cited to support the TI scenario of single-spin in-gap 
states, must have alternative explanations. Geometrical 
and final state effects are known to allow the detection of 
CD and spin-polarization in photoemission of non-chiral 
and non-magnetic solids [53, 54], and can prevent an un¬ 
ambiguous detection of spin-polarization asymmetries in 
YbBe, as discussed elsewhere [50]. 

Next we turn our attention to theoretical predictions 
of the YbBe electronic structure using the DFT method. 
We first reproduce the literature result [31] of a GGA ± 
SOC ± E (7 eV) calculation for YbBe in Fig. 3(a), which 
predicts a semimetallic band structure with a p-d band 
overlap at Ep. The local gapping at the band crossing 
points arises from rather weak bd SOC [55]. Since the p-d 
overlap anti-crossing points vary in energy around the X- 
point, the small local gapping cannot produce a full bulk 
gap, resulting in a complex semimetallic Fermi surface 
(FS), consisting of “lens” hole and “napkin ring” electron 
sheets. The calculated YbBe 4/ BE of 0.7-0.8 eV relative 
to the valence band maximum is in agreement with the 
experimental ARPES result in Eig. 1 of 1.05 eV which 
includes the 0.32 eV band gap. The location of the 4/ 
state far from Ep results in only a minor influence on the 
semimetallic ES that is thus very similar to predictions 
of the non-rare-earth divalent hexaborides [56-58]. 

In Eig. 3(b), we present an mBJ ± SOC ± E (7 eV) 
band result, overlaid with open-core pseudopotential sin¬ 
gle pass GW band result (dots) [34]. In both cases, the 
small p-d overlap of the GGA ± SOC ± U calculation in 
Eig. 3(a) is transformed into a small ^ 0.1 eV semicon¬ 
ductor gap with good quantitative agreement between 
the two methods [34]. This result clearly indicates that 
YbBe is a topologically trivial small band-gap semicon- 
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FIG. 4: (Color Online) (a) Pressure-dependent Yb valence 
state. For comparison, the experimental values are extracted 
from Ref. [33]. The mBJ + SOC U {7 eV) band structures 
(b) under P = 15 GPa, (c) under P = 30 GPa, and (d) under 
P = 50 GPa. 


ductor. Not surprisingly, slab calculations for both the 
YbBe (001) and (110) surfaces also show no topological 
in-gap surface states (see the Supplement [34]). 

This semiconductor result is reminiscent of the case 
of CaBe, whose early DFT-based semimetallic model for 
anomalous transport was revised to be that of a 1 eV 
semiconductor with the assistance of GW theory [59], and 
subsequently confirmed with ARPES [60] and other ex¬ 
periments using high-purity boron samples [61, 62]. This 
straightforward theoretical prediction for YbBe of being 
a topologically trivial semiconductor is in contrast to two 
recent calculations that predict YbBe to be a TI based 
on f-d band inversion [22] or p-d inversion [27, 30]. The 
flaws in these previous band calculations, resulting in in¬ 
correct Yb 4/ binding energies and mixed valency, are 
discussed in detail in the Supplement [34], along with 
angle-integrated valence band spectra from the ( 110 ) sur¬ 
face that provide definitive proof of the pure Yb divalency 
in YbBe [34], and thus additionally rule out these erro¬ 
neous theory calculations. 

A recent pressure dependent study of YbBe [33] ob¬ 
serves key results of (i) no structural transition up to 
30 GPa from x-ray diffraction, (ii) a rapid order-of- 
magnitude decrease in the resistivity up to 5 GPa, (iii) 
a pressure region of rather constant resistivity and Hall 
coefficient from 5 ^ 15 GPa, and (iv) a reemergence of 
thermally-activated resistivity above 15 GPa accompa¬ 
nied by a small increase in Yb valency from pure diva¬ 
lency to 2.09+. 

The theoretical calculation at 15 GPa in Fig. 4(b) 
shows a p-d overlap band structure and indicates that 
YbBe undergoes a semiconductor to semimetallic phase 
transition at an intermediate pressure. This occurs due 
to increase of p and d band widths and their wave func¬ 


tion overlap. Such a p-d gap to p-d overlap transition 
naturally explains the rapid initial decrease in resistivity 
with pressure, also observed in early pressure-dependent 
transport of YbBe [63]. A semimetallic state in the inter¬ 
mediate 5 ^ 15 GPa pressure regime is also suggested by 
the nearly constant Hall coefficient, which is attributable 
to a balance between electron and hole carriers [33]. This 
transformation to semimetallic behavior under pressure 
provides a further confirmation of the existence of a semi¬ 
conductor gap at ambient pressure where the ARPES 
experiments are performed. 

The theoretical electronic structures for even greater 
pressures of 30 GPa and 50 GPa in Eigs. 4(c) and (d) 
show an increasing p-d overlap such that the Yb Ad band 
ultimately touches the Yb 4/ band which remains at 
nearly the same BE. The Yb 4/ band exhibits only a 
small increase in bandwidth and slight centroid shift to 
lower BE but still remaining at the BE larger than 0.5 
eV. Nevertheless there is an increased mixing of Yb 4/ 
character into the p-states, as evidenced by the increasing 
band anti-crossing gapping that results from the Yb 4/ 
SOC interaction. The increasing Yb 4/ character above 
Ey implies a decreased /-occupation and mixed-valency. 
Quantitative analysis of the Yb valence under pressure 
is plotted in Eig. 4(a). The resulting mixed-valence, less 
than 10 % at the highest pressure, compares favorably to 
the experimental results derived from Yb I /3 x-ray ab¬ 
sorption measurements [33]. The experimental reemer¬ 
gence of a thermally activated resistivity (dRjdT < 0) 
above 15 GPa is plausibly due to the increasing 4/ SOC- 
induced local gapping, whereas the overall resistivity rise 
due to gapping is weakened due to the competition of the 
increasing p-d overlap and hence increasing hole and elec¬ 
tron ES volumes. The residual semimetallic conductivity 
can also explain the observed experimental low temper¬ 
ature resistivity plateaus [33]. 

This theoretical investigation allows us to comment 
generally on the feasibility of forming a TKI in actual 
materials. Since p-d states of opposite parity have inher¬ 
ently weak or negligible hybridization, the topologically 
non-trivial band inversion will have difficulty in forming 
a full insulator gap via hybridization alone. Therefore 
some additional external influence is required to open up 
an insulating gap of sufficient size to practically realize 
in-gap topological surface states. Here for the example 
of YbBe under pressure, the external infiuence is the hy¬ 
bridization mixing of the Yb 4/ states with the p-states 
and its larger 4/ SOC-induced gapping. However this 
effect is still too small for YbBe to develop a full BZ p-d 
overlap gap at experimentally achievable pressures. 

In conclusion, the fiat-band conditions of the non-polar 
( 110 ) surface allow ARPES measurements to definitively 
show that YbBe is a non-Kondo non-TI semiconductor, 
and it opens up a new method for the quantitative char¬ 
acterization of the bulk gap of other divalent hexaborides. 
This result is in good agreement with predictions of theo- 
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retical DFT+[/ calculations with proper treatment of 4/ 
correlations and inclusion of well-established gap correc¬ 
tion physics. Only under pressure does the topologically 
non-trivial p-d band inversion occur, but the system still 
retains a semimetallic electronic structure even up to high 
pressure beyond the onset of small Yb mixed valency. 
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Supplemental Material: 

Electronic structure of YbBs : is 
it a Topological Insulator or not? 


ARPES OF NON-POL AR (110) SURFACE 
A. Yb divalency and wide EDC 

The deep binding energy (BE) of the bulk Yb 4/ levels 
and lack of Yb^+ spectral weight near the Fermi level 
Ey is strong experimental evidence for pure divalency 
in YbBe- However, another important test for mixed- 
valency is the observation and quantification of the rel¬ 
ative amount of Yb^+ spectral weight in valence band 
photoemission measurements. Previous XPS and UPS 
measurements [1, 2] observed distinct Yb^+ contributions 
to Yb 4p core-level and valence band spectra, but the 
origin was ascribed to surface contributions. Similarly, 
a weak ferromagnetic moment corresponding to 2% bulk 
trivalency was initially reported in YbBe [3]. But later 
measurements showed stoichiometric YbBe to be diamag¬ 
netic, with surface ferromagnetism or paramagnetism ap¬ 
pearing only for non-stoichiometric samples [4]. 

In Fig. SI (a), we demonstrate the pure divalency of Yb 
from a wide valence band spectrum for a prepared YbBe 
(110) surface measured at hv = 120 eV. The spectrum 
exhibits Yb^+ 4 / 7/2 and 4 / 5/2 peaks that dominate over 
the very weak 4 - 12 eV BE range where Yb^+ spectral 
weight would be found. Overlapping in this energy range 
are B-sp valence band states whose relative cross section 
is low at this photon energy and which can be further 
suppressed using 5-polarization of the incident photons, 
as shown in the Fig. SI (a) inset. 

The line shape of the valence band structure, shown 
in the expanded view of this 4 ^ 12 eV range, exhibits 
no evidence for any Yb^+ spectral weight whose distinct 
/^^ ^ /^^ final state multiplet peak structure is illus¬ 
trated in the inset of Fig. SI (a) (dashed line). This 
example of a Yb^+ spectrum comes from a cleaved (001) 
surface that initially looks like the pure divalent (110) 
surface spectra, but then develops within only 4 hours 
the appearance of this Yb^+ signature, indicating the 
trivalent conversion of surface Yb atoms due to resid¬ 
ual gas adsorption. In contrast, the (110) surface with 
coplanar Yb and Be termination and 1x1 surface or¬ 
der, shown in Fig. Sl(b), exhibits no spectral changes 
after greater than 12 hours of measurement. Even if one 
integrates the B-sp spectral weight of the 5-polarization 
spectrum, after employing the background subtraction 
shown in Fig. SI (a), and assigns it to Yb^+, one obtains 
an upper bound contribution of less than 1%. This tight 
experimental upper bound on the pure divalency of Yb in 
YbBe immediately rules out both of the theoretical pre¬ 
dictions of f-d [5] and p-d [6, 7] band inversion scenarios. 
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FIG. SI: (a) Angle-integrated 120 eV spectra of YbBe (HO) surface demonstrating the pure divalency of Yb, and ruling out 
the calculations in Fig. S4(b,c). (Inset) Magnified view of the overlapping B-sp and Yb^^ regions showing the weaker B-sp 
cross section for s-polarization and the background subtraction for quantitative valence analysis. Also shown (dashed line) for 
comparison is the signature spectral line shape for Yb^^ ^ multiplet peaks acquired from an aged (001) cleaved surface 
of YbBe. (b) Low-energy electron diffraction of the prepared (110) facet showing 1x1 surface order, (c) Wide-energy 120 eV 
X-point spectrum of the YbBe (HO) surface showing the divalent Yb 4 / 7/2 peak at 1.05 eV and overplotted DFT bands along 
X-M. 


both of which predict mixed Yb valency. 

A 12 eV wide angle-resolved spectrum measured at the 
120 eV bulk X-point, shown in Fig. Sl(c), also illustrates 
the wide B-sp band dispersions in the 4 - 12 eV range, 
and the absence of any /^-independent Yb^+ peaks. The 
spectrum additionally illustrates the 1 eV and 2.3 eV 
binding energies of the Yb 4 / 7/2 and Yb 4 / 5/2 states, re¬ 
spectively. The Yb 4/ peaks are observed to be uniform 
across the prepared (110) sample surface, in contrast to 
the reported energy shifts between cleaves, and spatial 
inhomogeneity of the cleaved (001) surface [8-10]. The 
boron hole band that forms the valence band maximum 
is observed to disperse from below the Yb 4 / 5/2 state 
and significantly hybridize with the /-states as it passes 
through towards A small band width expansion of 
1.05 of the DFT bands was required to best match ex¬ 
periment. 


text, with an energy cut 0.3 eV below probing the 
valence band hole dispersion, shows two strong hole-band 
features ai ky = 0 for hu = 76 eV and 120 eV as well 
as at four weaker intensity second BZ X-points at ky = 

hiy = 120 eV was chosen for the kx — ky angular- 
dependent map in Fig. S2(b) with the same sample ori¬ 
entation. Similar to the photon-dependent map, strong 
valence band hole intensities are observed at three differ¬ 
ent X-points at ky = 0 and also additionally at four more 
X-points in the second BZ at ky = ±1.5A The bulk 
X-point at kx = —0.53A ^ was used for Fig. 1(e) in the 
main text and was checked for consistency at other bulk 
X-points. 


B. X-point location 

A suitable bulk X-point location for the band gap de¬ 
termination was chosen based on photon-dependent and 
angular-dependent mapping of the momentum-space. 
First a normal-emission photon-dependent kx — kz map 
was acquired with the sample (110) axis aligned verti¬ 
cally with the spectrometer slit. An energy cut at —1.1 
eV, between the two 4/ states, highlights a strong feature 
at 61 eV, as shown in Fig. S2(a) that aligns to the bulk 
Brillouin zone (BZ) using an inner potential parameter 
of 14 eV. 

Then the sample is rotated in azimuth to align the 
sample (100) axis to the vertical spectrometer slit and an 
off-normal photon-dependent map was acquired at fixed 
kx = —O. 53 A as indicated in Fig. S2(a). The ky — 
kz map at fixed kx^ presented in Fig. 2(a) of the main 



FIG. S2: (a) Normal-emission photon-dependent map with 
the sample (110) direction aligned to vertical spectrometer 
slit. The identity of the F-point defines the inner potential 
(Vb = 14 eV). The vertical line indicates the off-normal plane 
of the photon-dependent maps in Fig. 2(a) of the main text, 
(b) Angle-dependent map at 120 eV, showing the /c-location 
of the X-point spectrum used for the band gap analysis Fig. 
1 in the main text. 































C. K-dose 


dispersion can be written simply as [13]: 


To further characterize the conduction band minimum 
energy, surface K-dosing was used to create a small re¬ 
type band bending until the electron pocket was suffi¬ 
ciently deep to observe its non-parabolic dispersion. For 
this experiment, an X-point spectrum was obtained at 
76 eV instead of 120 eV due to the strong d-electron 
(p-hole) band intensity for s (p) polarization, allowing a 
good quantitative analysis of the non-parabolic conduc¬ 
tion (valence) band dispersion and semiconductor gap 
evaluation (see subsection D). 

Fig. S3 (a) shows the X-point p-polarization spectrum 
before K-dosing. Fig. S3(b) shows a time/dosing de¬ 
pendent energy distribution curve (EDC) cut at X il¬ 
lustrating the energy shifts of the valence and conduc¬ 
tion bands as well as the Yb 4/ state. The final EDC 
spectrum, shown in Eig. S3(c), is a sum of p- and 5- 
polarizations due to the weakening of the p-hole band 
with K-dosing. Because a rigid shift of both valence and 
conduction bands is observed in Eig. S3(b) with no dis¬ 
cernible surface band gap narrowing, we can claim that 
this K-dosed analysis is reflective of a bulk gap of mag¬ 
nitude 0.32 eV, based on non-parabolic dispersion visual 
fits to the polarization-dependent data that best high¬ 
light either the hole valence band or electron conduction 
band dispersions [11]. 


D. Non-parabolic Dispersion 


A two-band k • p model for non-parabolic valence and 
conduction bands in the presence of small band gaps [12] 
was used to assist in the determination of the valence 
band maximum and conduction band minimum energies. 
The two band k-p correction to the free electron parabolic 



- 0.4 0.0 0.4 0 2 4 6 - 0.4 0.0 0.4 

kx(A ') K-dose time (min) 


FIG. S3: (a) X-point spectrum at hi> — 76 eV of the YbBe 
(110) surface before K-dosing showing Ep pinned to the bot¬ 
tom of the Yb bd conduction band, (b) K-dosing EDC image 
of the valence and conduction band energy shifts, (c) X-point 
spectrum after K-dosing allowing conduction band dispersion 
analysis for band gap determination. 


e{k) 


/ 4 

2 yY Eg 2m* 



( 1 ) 


where Eg and m* are the energy gap parameter and the 
electron effective mass, respectively. 

The energy gap dependence of Eq. (1) exhibits 
the well-known progression towards a purely linear V- 
shaped dispersion at zero gap, corresponding to a Dirac 
semimetal. The valence band dispersion in Eig. 1(d) of 
the main text exhibits linearity away from a distinctly 
rounded band maximum and is nicely fit using the nom¬ 
inal experimental gap value of 0.3 eV and an effective 
mass of m* = 0.12. This dispersion shape for interme¬ 
diate small gap values contrasts sharply with topological 
Dirac cone surface states in Bi 2 X 3 (X = Se, Te), where 
the dispersion is mostly linear near the Dirac point and 
deviates further away (e.^., hexagonal warping). 

The conduction band dispersion, revealed after K- 
dosing in Eig. S3(c), is similarly fit to Eg = 0.3 eV and 
m* = 0.1. 


COMPUTATIONAL DETAILS 

The band calculations were performed by us¬ 
ing the full-potential linearized augmented plane- 
wave (ELAPW) band method, as implemented in 
the WIEN2K package [14]. Eor the DET calcula¬ 
tions, the PBE (Perdew-Burke-Ernzerhof) exchange- 
correlation functional was used in the GGA (generalized- 
gradient approximation). The spin-orbit coupling (SOC) 
was included in a way of the second variational method. 
In the GGA + SOC + U method, a correlation energy of 
U = 7 eV was chosen to obtain the correct experimental 
value of the Yb 4/ BE of ~ 1 eV, which is consistent with 
the previous calculations [15, 16]. The mBJ potential cor¬ 
responds to an orbital-independent semilocal exchange 
potential that mimics the behavior of orbital-dependent 
potentials and has been shown to provide band gap cor¬ 
rections in good agreement with the improved many- 
body but more computation-demanding GW calculation 
[17, 18]. A single pass open-core pseudopotential GW 
calculation was also performed using VASP [19], which 
effectively removes the Yb 4/ levels from the calculation. 
This is reasonable in view of the deep 4/ BE in YbBe- 

In the GGA U {U: on-site Coulomb interaction 
of Yb 4/ electrons) calculations, the scheme of a fully- 
localized limit was used for the 4/ electrons. The BZ 
integration was done with al7 x 17 x 17 /c-mesh, and, 
in the ELAPW calculations, muffin-tin radii of Yb 

and B were chosen to be 2.50 and 1.59 a.u., respectively, 
and the plane-wave cutoff was RMT^max = 7. 

Eor the GW calculation, a single pass open-core 
pseudopotential Go Wo calculation was performed using 
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GGA + SOC + U (7 eV) GGA + SOC + U (4 eV) GGA + SOC + U (0 eV) 



FIG. S4: DFT-Wien2k band structures of YbBe. (a) GGA + 
SOG + (7 eV) calculation yields a semimetallic p-d overlap 

with anticrossing gaps, (b) GGA + SOG + (4 eV) calcula¬ 

tion yields SOG enhanced gapping of overlapping p-d states, 
similar to the VASP calculation by Neupane et al. [6]. (c) 
GGA + SOG + f/ (0 eV) calculation yields f-d band inver¬ 
sion, similar to the LDA + Gutzwiller calculation by Weng et 
al. [5]. 

VASP with a 4 X 4 X 4 /c-mesh, 400 unoccupied bands 
and 80 mesh size in the omega frequency. 

YbBe crystallizes in the CsCl-type structure, like 
SmBe. The ambient pressure calculations were per¬ 
formed using the experimental YbBe lattice parameter 
of a = 4.1792 A and internal crystal parameter xb = 
0.202 that specifies the relative inter-octahedron B-B 
bond length [15]. For pressure-dependent calculations, 
a linear relation between the lattice constant and ap¬ 
plied pressure up to 30 GPa was established using x-ray 
diffraction data in Ref. [20], before starting the xb re¬ 
laxation. To obtain the structure information under 50 
GPa, where the lattice constant was not accessible in 
Ref. [20], the fully relaxed calculation was performed. 
All electronic structure calculations under pressure were 
done using the mBJ + SOC + 7/(7 eV) scheme. We 
have checked the pressure dependent U value by a con¬ 
strained LDA calculation [21] and found that U ^ 7 eV 
is obtained over the whole pressure range, which could 
be interpreted as showing the strong localization of the 
Yb 4/ orbital. 

DFT + SOC + U BAND STRUCTURES FOR 
YbBe WITH VARYING U 

Two recent calculations predict YbBe to be a topo¬ 
logical insulator (TI) based on f-d band inversion [5] or 
p-d inversion [6, 7]. The first was an LDA + Gutzwiller 
calculation with U = 6 eV that erroneously predicted a 
topological f-d inverted band structure very similar to 
SmBe [22], with a highly incorrect 4/ BE of ^ 0.1 eV. 


The LDA + Gutzwiller result [5] is empirically very sim¬ 
ilar to a 7/ = 0 DFT result in Fig. S4(c) in terms of the 
4/ BE and the f-d band inversion at Ep. 

A second proposed topological band structure for 
YbBe, based on DET + U pseudopotential VASP cal¬ 
culations using small U values of 4 eV [6] or 5 eV [7], 
obtains a p-d band overlap but with a small insulating 
gap throughout the bulk BZ, similar to that in Eig. S4(b) 
[23]. The effect of the smaller U values is a smaller Yb 4/ 
BE, 0.3 (0.5) eV for 7/ = 4 (5) eV, whose closer proxim¬ 
ity to the p-d overlap region at Ep results in significant 
hybridization of the Yb 4/ and B 2p states. This allows 
the Yb 4/ SOC to open a larger p-d gap that otherwise 
is intrinsically tiny due to the smallness of the bd SOC. 
The hybridization of Yb 4/ states into the p-d overlap 
region also results in a mixed-valent Yb 4/-occupation 
evaluated to be n/ ^ 13.8, which does not agree with the 
experimental divalent state of Yb. Chang et al [7] also 
reported that the SOC-induced gapping of the p-d over¬ 
lap at 7/ = 4^5 eV evolves into a topologically trivial 
p-d gap above U = 6.5 eV. This result is in contradiction 
to the earlier literature result [15] and our Wien2k result 
shown in Eig. S4(a). As shown in Eig. S5, the DET + 
U calculations for increasing U yield only a small grad¬ 
ual decrease in the p-d overlap energy without opening 
up a semiconducting p-d gap even for unphysically large 
U values as high as 15 eV. This is consistent with a de¬ 
creasing influence of the Yb 4/ states on the p-d states 
as they are pushed farther away from Ep. Namely, large 
U is not the formative parameter for the semiconductor 
gap. 

One might consider that the U > 6.5 eV VASP predic¬ 
tion of a topologically trivial semiconductor is in agree¬ 
ment with the Wien2k mBJ + SOC + 7/ semiconductor 
result in Eig. 3(b) of main text. However this is en¬ 
tirely fortuitous because the VASP calculation does not 
contain the band gap renormalization treatment of the 
mBJ or GW calculations in Eig. 3(b). Chang et al. [7] 
also presented a GW calculation for YbBe but for U = 
0 with the goal of demonstrating that GW alone cannot 
produced the p-d band inversion. In contrast, the VASP 
GW calculation in Eig. 3(b) is an open-core calculation, 
which simulates a large effective U that removes the Yb 
4/ states far from E’p. Thus the GW calculation in Eig. 
3(b) treats the p and d correlations that renormalize the 
p-d overlap into a small semiconductor gap. 

We have also checked DET + U electronic structures 
using VASP. The crystal structure was fully relaxed us¬ 
ing the GGA + SOC scheme (a = 4.13 A, = 0.201). 
The resultant crystal structure was quite similar to the 
experimental structure. Band structures of YbBe were 
calculated with various U values in both relaxed and ex¬ 
perimental structures, and it was found that the band 
structures for the two cases are almost the same. We have 
found that our VASP results are quite consistent with our 
Wien2k results. There is a sensitivity of the p-d overlap 
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(b) 



U (eV) 


FIG. S5: (a) The GGA + SOG + U calculations using Wien2k, showing the persistence of the p-d semimetal overlap up to 
unphysically high values of U. Red, blue, and green color symbols represent Yb 4/, Yb 5d, and B 2p characters, respectively, 
(b) Summary comparison of the p-d band overlap and the binding energy of Yb 4/ band (^ 4 /) between our Wien2k and the 
VASP calculations by Ghang et al. [6, 7]. The SOG-induced p-d band gap by Ghang et al. (SO gap) are also plotted. 


to the internal structural parameter xb [24]. One could 
speculate that differences in structural input could be a 
possible source of the discrepancy between our Wien2k 
(and VASP) and the published VASP calculations [ 6 , 7], 
but no structural information (lattice constant or xb) 
was provided in Refs. [ 6 , 7]. 

Figure S5(b) shows the comparison between the VASP 
results by Chang et al [7] vs. our Wien 2 k results, which 
highlights the discrepancy for 7/ >6 eV. We can note 
discrepancies between two in (i) the binding energy of Yb 
4/ band (^ 4 /) and (ii) p-d overlap. The VASP results 
of the latter by Chang et al. look merely shifted Wien 2 k 
results. 

In summary, we note that the incorrect f-d overlap [5], 
and p-d overlap [ 6 , 7] calculations did not cite the earlier 
Wien 2 k calculation of Ref. [15], and hence failed to com¬ 



FIG. S6: DFT-VASP open-core slab calculations for YbBe 
(001) surface. For the purpose of getting a trivial p-d gap 
without the mBJ scheme, a lattice constant of YbBe is en¬ 
larged to be a = 4.5A. (a) Be termination, (b) Yb termina¬ 
tion. Both of terminations give topologically trivial surface 
states. The surface states (red lines) below Ep are mostly of 
B-p character, while those above Ep in (b) have mixed B p 
and Yb d characters. 


pare and discuss their discrepancies with the published 
literature. We also note that the complicated intermixing 
of all three p, d and / states near Ep with interplay of 4/ 
SOC and mixed valency discussed in the main text for 
the high-pressure phase of YbBe is the same physics at 
the heart of the erroneous claim of a TI p-d overlap gap 
at ambient pressure [ 6 , 7], albeit errantly realized from a 
too small value of U. 

A remaining open question about the ambient pressure 
phase is the quantitative disagreement between the ex¬ 
perimental 0.3 eV and theoretical 0.1 eV gaps, which 
would be worthwhile to be investigated theoretically. 
Further theoretical investigation of the high pressure 
phase is warranted to address whether any low energy 
scale dynamical correlations {e.g. Kondo effect) can 
emerge from the small mixed valency, or whether any 
non-trivial topological p-d parity inversion effects can be 
realized while in a partially gapped semimetallic state. 


DFT SLAB CALCULATION FOR YbBe (001) 
SURFACE 

In order to examine the surface states, we have also car¬ 
ried out DFT open-core slab calculations for the YbBe 
(001) surface using VASP. We have used slab structures 
composed of 9 Yb and 8 Be layers for Yb termination 
and 8 Yb and 9 Be layers for Be termination, and for 
both terminations, ^ 30 A vacuum layers were consid¬ 
ered. We have first relaxed the slab structures for both 
terminations, and then self-consistent calculations were 
performed with a 10 x 10 x 1 /c-mesh. For the purpose 
of getting a topologically trivial p-d bulk gap in the DFT 
+ SOC + 7/(7 eV) without employing mBJ, we em¬ 
ployed an increased lattice constant of a = 4.5A in the 
slab calculations. 

As shown in Fig. S6, surface states (red lines) ap¬ 
pear for both Be and Yb terminations. These sur- 
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face states, however, are just topologically trivial surface 
states, which are often observed in the hexaboride sys¬ 
tems [25, 26]. The Dirac-cone-like surface states observed 
in ARPES are not obtained. In fact, such Dirac cones 
are expected to be quantum well states coming from the 
band bendings in the polar (001) surfaces with different 
charge terminations [10, 11]. 

The DFT open-core slab calculations for the non-polar 
(110) surface also produce similar surface states, one deep 
surface band of B p character and two surface bands near 
Ey of mixed B p and Yb d characters that arise from 
broken bond states from two surface B atoms [27]. All 
of those surface states are topologically trivial. 

LINK BETWEEN YbBe AND SmBe 

There are strong differences in issues of the two mate¬ 
rials YbBe and SmBe ie., pure-divalency and p-d over¬ 
lap/semiconductor gap issues in YbBe and mixed valency 
and f-d hybridization gap issues in SmBe. Also note that 
while YbBe (001) ARPES shows “false” in-gap states 
that are actually band-bending quantum well states, the 
YbBe (110) ARPES show no states within the 0.3 eV 
semiconductor gap, thereby disproving the TI scenario 
for YbBe. 

In contrast for SmBe the 20 meV hybridization gap 
exhibits similar ARPES in-gap states for both the polar 
(001) and non-polar (110) surfaces [28]. The existence 
of the X-point in-gap states for SmBe is insensitive to 
the different polar surface terminations, and similar in¬ 
gap states are observed for the non-polar (110) surface 
too [28]. Hence the polar nature alone cannot explain 
the in-gap states, although polar effects are still being 
discussed for the (001) surface, e.g., in the context of 
2D-hke d-states outside the gap [29]. 

We have also checked the mBJ scheme for SmBe. The 
reason why we used the mBJ scheme for YbBe is to con¬ 
sider the correlation effects of delocalized Yb 5d and B 2p 
electrons near Ep. While it is very crucial for YbBe that 
the mBJ correction changes the p-d overlap into a p-d 
trivial gap, for SmBe, the Sm 4/ bands are dominating 
near Ep, and so the mBJ method does not play a crucial 
role. In fact, we have confirmed for SmBe that the DFT 
and mBJ results are essentially the same [27]. 

For SmBe, dynamical mean-field theory (DMFT) is es¬ 
sential to address the dynamic energy correlations that 
result in strong low energy scale renormalization of the /- 
states near Ep as well as the mixed valency. In contrast, 
for YbBe, with the large binding energy of the 4/ states 
and pure divalency, one would not expect any low en¬ 
ergy scale dynamical correlations to emerge from DMFT 


calculations. Indeed we have checked and confirmed for 
YbBe that DFT + DMFT gives essentially identical re¬ 
sults as DFT + 1/, be., a proper 4/ binding energy and 
an uncorrected p-d overlap [27]. The latter is because the 
correlation effects of delocalized p and d electrons are not 
treated within DMFT. 
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